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On Brownian motion on the plane with 
membranes on rays with a common endpoint 
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o 

' Abstract 

04 ■ We consider a Brownian motion on the plane with semipermeable 

^ ' membranes on n rays that have a common endpoint in the origin. We 

I obtain the necessary and sufficient conditions for the process to reach the 

. origin and we show that the probabihty of hitting the origin is equal to 

_ ' zero or one. 

! Introduction 

^ . Let Ci , . . . , Cn he given rays in such that they share a common 

! endpoint. Suppose polar coordinates in are denoted by (r, ^p),r > 
0,ipe [0, 27r) and 



Ck = {{r,'f) : r > 0, (f = cpk}, 



where < v^fc < 2TT,(pj ^ i^k,3 k,j,k E {l,...,n}. 
^ I Let rifc, k = 1, . . . ,n, denote the unit vector normal to that 

O ■ points anticlockwise and let Vk be a vector in such that {v^, rtk) = 

1. The angle between and Vk denoted by 6^^ G (— |, |) is referred 
^ ! to as a positive if and only if points towards the origin. The case 

Q\ [ of 77, = 5, 6'i > 0, 6*2 > is shown in Figure [TJ 

O ■ Let 7fc, /c = 1, . . . , n, be a real constant such that |7fc| < 1. 

Consider the stochastic differential equation in 

X 



dx{t) = dw{t) + -fkVkdLl" (t) (1) 



k=l 



with initial condition a;(0) = x^, x° G M^. 

In this equation {w(t))t>o is a Wiener process in M?, (L^'=(t))t>o 
is a local time of the process {x(t))t>o at the ray Cfc, and it is defined 
by the formula 

Ll-{t) =\im^ I (2) 
eio 2e Jq 



^Research partially supported by Ministry of Education and Science of Ukraine, Grant 
F26/433-2008, and National Academy of Sciences of Ukraine, Grant f04-2008. 
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Figure 1. 



where = {x E M.'^ : 3y E Ck, s G [—1, 1] such that x = y + esuk}- 
The solution of this equation can be regarded as a Wiener process 
in skewing on the rays Ci, . . . , c„ until the time it reaches the 
origin (see Section 1). The skew on is defined by 7^ being a 
coefficient of permeability. When 7^ is equal to 1 we have a Wiener 
process reffecting instantaneously at Ck in the direction of Vk- If 7fc 
is equal to —1 the process reffects at Ck in the direction of —Vk- If 
7fc G (—1, 1) a semipermeable membrane is on c^. 

The aim of the investigation is to show that the process hits the 
origin, when starting away from it, with probability zero or one. We 
obtain an explicit expression in variables 9i, . . . , 9n, 71, ... , 7n, ^1, . . . , ^, 
whose sign this probability depends on. 

The case of the wedge (two rays) with reffection on its sides 
was studied by Varadhan and Williams [Tj. They proved (in our 
notation) that the process does not reach the origin if 6'i — ^2 < 
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and does reach the origin if — ^2 > 0. In this case our result 
coincides with that of Varadhan and Wilhams (see Example 2). 



1. Construction of the process 

In this Section we build the process which is a solution of equation 
([1]) up to the time it hits the origin. 

Conventionally, we put Lfni+k = ^k,Cni+k = Ck, jni+k = Ik, k = 
l,...,n, I e N. Define ^k = <^k+i - Vfc- Let Sk = {{r,ip) : r > 
0, ^Pk ^ ^ V^fe+i}, A; = 1, . . . , n, be the wedge CkOck+i- Let (pi = 0. 

Remark 1. Without loss of generality we can assume that C,k < 
n, k = 1, . . . ,n. Indeed, otherwise we can introduce additional ray 

Cad = {{r,^) : r > 0, (/^ = tt}, 

and put 7ad = 0. 

So, throughout this Section we assume < 7r/2, A; = 1, . . . , n. 

Lemma 1. Given x*^ G M^, there exists a unique strong solution of 
the equation (QP up to the first time of hitting the origin. 

Proof, (i) The process (a;(t))f>o is an ordinary Wiener process up to 
the first time of hitting one of the rays ci, . . . , . Define 

Ti = inf {t > : x(t) G Cj^ for some ji G {1, . . . , n}} . 

Without loss of generality, let x{ti) G Ci. Until the process reaches 
another ray only the local time at the ray Ci can increase. Put 

T2 = inf {t > n : x{t) G for some js G {1, . . . , n}, ja 7^ ji} • 
Then for all t < T2 equation ([1]) has the form 



where {xi{t),X2(t)) are the cartesian coordinates of the process {x(t))t>o, 
{wi(t))t>o and {w2(t))t>o are independent one-dimensional Wiener 
processes. Note that by Remark [1] 



dxi(t) 
dx2{t) 



dwi it) — 7i tan OidL^^ (t) , 
dw2{t)+-fidLl^{t), 



(3) 
(4) 




'S.[-e,e]{x2{s))ds, t < Tg, 
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where {L'^^{t))t>o is a usual local time of the process (x2(t))t>o at 
zero. There exists a unique strong solution of equation (jlj) (cf. 
Substitution of solution to (j4]) into ([3]) allows us to get the solution 
of the latter. Hence we have proved existence and uniqueness of 
solution to equation ([T]) up to the time T2. 

(ii) Put for m > 2 

Tm = inf {t > Tm-1 : x(t) e Cj^ for some jm G {1, . . . , n}, 7^ jm-i} 

Assume x{t2) G C2. Then for T2 < t < T3, the equation ([1]) has the 
form 

dx,{t) = dw,{t)-^2^-^^^^^^±^dLl-{t), (5) 

cos6'2 

dx2{t) = dw2it)+^2^^^^^^-i^dL:'it). (6) 

COS U2 

We do the clockwise rotation of coordinates defined by the formulas 

x[ = xi cos^i + X2 sin^i, 
x'2 = — xi sin,^i + X2 cos,^i. 

This map takes the ray C2 to the ray ci, i.e. C2 = ci. Moreover, it 
preserves the angles. Using the Ito formula we get the stochastic 
differential equations 

dx[{t) = dw[{t) --f2tane2dL''},{t), (8) 
dx'2{t) = rf^(t) +72rfLj(t), (9) 

which holds true for T2 < t < r^. Here w[(t) = cos^idwi(t) + 
sm^idw2(t), t>0, and W2(t) = — sin ^idwi(t) + cos ^idw2{t) , t>0, 
are independent one-dimensional Wiener processes as a square inte- 
grable martingales with characteristics t. By the arguments similar 
to those of (i) there exists a unique strong solution to equations ([H]) 
and ([2]) for r2 < t < It means, because the inverse transformation 
exists, that the equations ([5]) and ([6]) has a unique strong solution 
for T2 < t < T3. The case of x{t2) G c„ can be treated analogously. 
Thus we have built a strong solution to equation ([1]) up to the time 
Ts and proved the uniqueness of it. 

(iii) Repeating this procedure of localization of equation ([1]) we 
build a solution to ([1]) up to the time ( = limm^ooTm- It is easily 
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seen that the process (x{t))t>o does not reach infinity in finite time 
a.s. and does not intersect any ring a < \x\ < b, < a < b,x ^ M^, 
infinitely many times on [0, T]. So, ( = +00 or x{(—) = 0. □ 



2. Coordinate transformation and change of time 

As above, we assume that < 7r/2, k = 1, . . . ,n. 

To investigate the value of probability Pa;o{C < 00} we make 
certain transformations of the process. 

Let us change coordinates by the formulas 



where 



Therefore, 



Xl = — In yxf + X2 = — Inr, 

r cos ip = Xl, 
r sin ip = X2- 



-e^^cosx2 — e^^sinx2 
-e^^ sin X2 e^^ cos X2 



Xl = e C0SX2, (10) 
X2 = e~'^^sin5;2. (11) 

The Jacobian matrix of this transformation is 

^(^1,^2) _ / d{xi,X2) \ _ 

d{xi,X2) \d{xi,X2)J 

Moreover, Axi = Ax2 = 0. 

The image of the ray Ck,k = l,...,n, under this map is the 
straight-line 

Ck = {(^1,^2) : Xl e M,X2 = ipk}- 
Clearly, the image of 5*^ is the strip Sk = {{xi,X2) : Xi e M, < 
X2 < <Pk+i}- As above we put ipni+k = ^k, Cni+k = 4, k = 1, . . . ,n,l e 
N. 

Note that if the function / is twice continuously differentiable on 
\ {0}, then by the Ito formula 

n 

df{x{t)) = {Vf,dw{t))+J2^k{Vf,Vk)dLl''{t)+-Af{x{t))dt, t < (. 

k=i 

(12) 



5 



So 

dxi{t) = -e^^^^\cosx2{t)dwi{t) + sm.X2{t)dw2{t))- 
'^ (13) 
_e-iW (t,^, cosx2(t) + i;fc2 sinSaW) c^i^^H^), ^ < C, 

k=l 



where {vki,Vk2) are the cartesian coordinates of the vector Vk- As 

easily seen, Vki = - """cos ffc2 = ^^^fe^^- '^^^ process {Ll^{t))t>o 
increases only at those moments of time when x E and, conse- 
quently, X2 = fk- Besides, the process Wi(t) = — cosx2(t)dwi{t) — 
smx2(t)dw2{t),t > 0, is a Wiener process as a square integrable 
martingale with its characteristic being equal to t. Thus, we can 
rewrite the equation (fT3|l as follows 

n 

(ixi(t) = e*i(*)(iwi(t) + e^i(*) ^7fctan4c;L^'=(t). (14) 

fe=i 

By analogy, 

n 

dx2{t) = e*iWc/w)2(t) + e^'^(*) ^7fcrfL^H^), (15) 

k=l 

where W2it) = cosx2(t)dw2it) — smx2(t)dwi{t), t > 0, is a one- 
dimensional Wiener process independent of the process (wi(t))t>o- 
From f[T^ . IHM we get the equation 

n 

dx{t) = e^'^'^dw{t) + e^^W ^7fcCfcdL^^(t), (16) 

k=l 

where w(t) = {wi(t),W2(t)),t >0;Vk = (tan6'fc, 1), = 1, . . . , n. 

We emphasize that L^,* is not a local time of the process x but 
that of the process x. Let us give a form of equation f|T6l) that 
involves L^*" being a local time of the process x at the line Ck instead 
of L^*^, /c = 1, . . . , n. 

To do this we make use of the following lemma. 

Lemma 2. Let u{t), t G [0,T], he a continuous function on [0,T] 
such that for all t G [0, T] there exists a limit 

1 /■* 
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Suppose f is a continuous function on [0, T] with positive values, 
g{s,e) is a function on [0,T] x [0, T] with non-negative values such 
that g{s,e)/e — > /(s) as e I uniformly in s. Then for all t G [O.T] 

^^^T' I ^[-~9{s,e),g(s,e)\{u{s))ds = [ f{s)dL{s). 

Proof, (i) The statement of Lemma is easily justified for g{s,e) = 
el[[a^h]{s), where a,b are constant, < a < 6, and, consequently, it 
is true for linear combinations of such functions. 

(ii) Let fn, fn be sequences of piecewise constant functions on 

[0,T] such that fn < f < fn and fn -> f, fn ^ f as n ^ oo 
uniformly. By (i) we have 

[ fn{s)dL{s) = \im^ [ l[-eMs),sMs)]{u{s))ds < 

<lim;^ / ll[^ef{s),e.f{s)]{u{s))ds <lim^ [ < 
£iO Jo Jo 

Passing to the limit as n — > oo we get 

^^'^'^ I ^[~ef{s),ef{s)]{u{s))ds = f f{s)dL{s). 

^■10 ^£ Jo Jo 

(iii) Since — f{s) as e | uniformly in s G [0,T], then for 
each t G [0, T], 6 > there exists es > such that for all e < es the 
inequalities 

e{l-6)f{s)<g{s,e)<e{l + 6)f{s) 
hold for all s G [0,t]. Similarly to (ii) we have 

(l-'J) / JXs)dL{s) <ljm^ [ ]I[_g(,,^),g(,,^)](M(s))rfs < 

Jo slO Jo 

-^^hJo '^l-9i^'^)Ms,e)]{u{s))ds < {1 + 6) f{s)dL{s). 
Letting 5 | we obtain the statement of the Lemma. □ 
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Now let us express L^* in term of L^*. 
Lemma 3. The process 

Jo 

is a local time of the process {x(t))t>o on the straight-line Ck- 
Proof. Let k = 1. According to ffTTl) 

L2 (t) =\im^^ I\[^,^4x2{s))ds = 

1 /"* 

= ^^f?^ / 2[-e-*iWsme,e-*i(»)sine]S{cos£2(.)>0}(e"'''^'^sinX2(s))rfs 

1 /•* 

= lim^ y Ihe-*lWsm£,e--lWsm£](a;2(s))cis = 



t 



The last equality follows from Lemma El This proves Lemma [3] in 
the case of /c = 1. The statement of Lemma for /c = 2, . . . , n, can 
be obtained analogously. □ 

Using this Lemma we get the form of equation (fT6!) involving 
only the process x as follows 

n 

dx{t) = e^^Wciw}(t) +e2^^W^7fc^^fcrf4'W- (17) 

k=l 

Let us make a change of time. Put A(t) = j^e^^^'^'^Ms, t < 
and denote rj = e'^^^^^^ds. Define the process 

X{t) = x{A-\t)), t < 7], 

where ^ 

A~\t)=mf{s : [ e^^'^'^^duyt}. 
Jo 
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The equation (fTTI) implies the stochastic integral equation 

x{A-\t)) = £(0) + / e*^«dil;(t) + 



71 „ 

k=l "^0 



Similarly to [S], Theorem 7.2, it can be proved that the process 
w(t) = e^^^^''dw{s) is a Wiener process in up to the time 

V- 

Lemma 4. The process 

^0 

is a local time of the process {x(t))t<n on the straight-line Ck- 
Proof. For k = 1 we have 



= liml !\^^,^,ii^{A-\s)))ds - 
= lim- / ]I[„,,,](x2(^))e^^^(^)c? 

elO IE Jo 



X 



The cases of /c = 2, . . . , n can be treated analogously. □ 

Thus using (ITSl) we get the stochastic differential equation for the 
process (x(t))t<^ as follows 

n 

dx{t) = dwit) + J2^kVkdLf{t). (19) 
fc=i 

Remark 2. The process (x(t))t<^ is a continuous Markov one as a 
result of the random change of time for the Markov process {x{t))t>o 
which is continuous (cf. [2]). So, either x{t) is defined for all t > or 
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limf^^_ \x{t)\ = oo. It is easy to see that any solution to SDE f|T9|) 
does not blow up in finite time a.s. Moreover, in the next Section 
we will show that EL^^(t) < oo for each t > (see fl22l) ). Then, by 
(fT9l) . the trajectories of x{t) do not reach infinity in finite time a.s. 
Thus the process x[t) is defined for all t > 0, i.e. P{?7 = 00} = 1. 

3. Ergodic behavior of some Markov chain re- 
lated to the Brownian motion with membranes. 

For A; G N, let pk = ^-^,qk = Recall that = V^fc+i 

. . Tm,m G N, is defined in the proof of Lemma [H Put ^0 '■= ^n- 

The process X2(t) behaves on (r™, Tm+i), m = 1, 2, as an or- 
dinary skew Brownian motion with skewing at the point X2(Tm) G 
{ipi, . . . ,(Pn} (Fig. E]). The skew Brownian motion is a strong 
Markov process as a continuous Feller process (cf. [2], Theorem 
3.10). This yields {x2{Tm))m>i form a homogeneous Markov chain 
with the phase space {(fi, . . . , (fn} and the transition matrix 

'0 pi .... qi' 

q2 P2 .. 



. . Qn-l Pn-1 
Pn ■ ■ Qn _ 

where 

Pk = P^fe {x2{(^^k-i = Vk+i] ,qk = 1-Pk, 

Oa = inf{t > : X2(t) = a}. 
One can see (cf. [3], Ch. 3) that, for /c = 1, . . . , n, 

_ _ Pk^k-i 

Pkik~i + qkik ' 
~ _ Qk^k 

Pkik~i + qkik ' 

Now we find M^. := E^j.L|^ (cr^j._j Aa^j.^J. Equation f|T9l) implies 

E^fe£2(c^^fe_i A a^.^J = <^fc + 7fcE^,L£^^( 



(20) 
(21) 



10 



Figure 2. 



Using ([20]),([2ll) we get 

Mfc = . (22) 

3.1. Assume < pk < 1? k = l,...,n. The Markov chain 
{x2{'Tm))m>i IS irreducible and all of its states are positive-recurrent. 
This implies the existence of a unique invariant distribution (vri, . . . , vr^ 
which is the solution to the following system of equations 

TTfc = Tik-ipk-i + TTk+iqk+i, k = 2, . . . ,n - 1, (23) 

TTl = Vr„p„ + 712^2, (24) 

n 

Y^Tik = 1. (25) 

k=l 

Let us find it out. We can rewrite the equation (1231) in the form 

T^kiPk + Qk) = TTfe-iPfc-i + nk+iQk+i, k = 2,...,n-l. 
This implies the equality 

TTfc+igfc+l - TTkPk = TTkQk - TTk-lPk-l 

valid for /c = 2, . . . , n — 1. Put 

C := TTk+iqk+i - T^kPk, k = 1, . . . ,n - 1. (26) 
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Let k = 1. Then from (l26l) we have 



vr2 = C^ + 7ri?i. (27) 

g2 q2 



Put /c = 2 in fl26l) . Taking into account fl27j) we obtain 

yiys g2g3/ q2q3 

Finally by recursion we arrive at the formula 

^, 1 , Pk-l . P2-...-pk-l\ , Pl-...-pk-l 
TTk = C { — + — H h — : — + TTi — 



gfc qk-iQk q2 - ■■■■qk J q2 - ■■■■qk 

justified for /c = 2, . . . , n. 

Substituting this expression for 7r„ in the equality 

VTigi - -KnPn = C, 

which is a simple consequence of (l24l) . we get 

1 + r„ + rn-i^„ H h r2 ■ . . . ■ f„ . 

71"! = — z C, ltpi-...-pnt qi-----qn, 

where fk = ^,k = l,2,...,n. Note that pi ■ . . . ■ p„ 7^ gi ■ . . . ■ g„ if 
and only if ■ . . . ■ 7^ gi ■ . . . ■ g„. Taking into account that we 
can chose the ray Ci in an arbitrary way, we get for k = 1,2, . . . ,n 

1 + r-k^i + fk^2fk~i H h ri ■ . . . ■ fk-ifk+i ■ ■■■■fn „ 

if Pi ■ . . . ■ p„ 7^ gi ■ . . . ■ g„, 

One can find the value of C by substitution of the last formula 
into (|25l) . Computations lead to the formula 

vTfc = CiPk, (28) 

where 



^j=0 ^fc+l ■ • • • ■ qk+jPk+j+l ■ ■ ■ ■ ■ Pn+k-l^k+1 ■ ■ ■ ■ ■ ^k+j^k+j " • • • " ^n+k- 

nr=l, i^fcfe"^*-! + 9i6) 

n 

fc=l 



2 
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li Pi Pn = qi ■■■■■ Qn, it is easy to see that 

1 

I -U El -L. P'^P2 I ... I Pl-.-.-Pn-l ' 

Pl-...-Pk-l , „ 

TTi— z — , K = 2, . . . ,n. 

q2 - ■■■■ Qk 

Now suppose < pfc < 1, k = 1, . . . ,n, oi < p^ < 1, k = 
1, . . . , n. It is easy to check that formula fl28p has no singularities in 
this case. 

Thus we have proved the following statement 

Lemma 5. Let, for k = < Pk < 1 or let, for k = 

1, . . . ,n, < pfc < 1. Then there exists a unique invariant distribu- 
tion {7ik)^^i of the Markov chain (x2(rm))m>i- 

If Pi ■ . . . ■ Pn ^ qi ■ ■ ■ ■ ■ qn this invariant distribution is given by 
formula [^) . k = 1, . . . ,n. 

If Pi ■ . . . ■ Pn = qi ■ ■ ■ ■ ■ q-n the invariant distribution is given by 
formula / f^) . /c = 1, . . . , n. 

Remark 3. The method of finding the invariant distribution was 
proposed by Professor Alexey Kulik. 

Now partial cases when some of pk are equal to and some of 
them are equal to 1 are left unattended. We consider them in Sub- 
sections 3.3, 3.4. 

3.2. Further on we make use of the following well-known result. 
Given a homogeneous Markov chain with state space {0,1,. . . , N} 
and transition matrix 

'O 1 

xi 1/1 ... 



... xn^i vn^i 
1 

where < < 1, < < 1, i = 1, . . . , - 1. 

As is known there exists a unique invariant distribution 

{fi]i=o = {f{xu • • • , xn-uVu • • • , Vn-u A^)}^o 



TTi = 
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where 



•^0 -[^ I 1 I ^/i I I yi-----yN-2 I yi-.-.-VN-i ' ('^^) 

Xl X1X2 XI-...-XM-1 X1-...-XJV-1 

/i = -/o, (31) 

/. = ^^^——^fo,^ = 2,...,N-l, (32) 
• . . . • 

/tV = VN-lfN-l- (33) 

3.3. Suppose pfco = 1, pfc^ = . . . = = 0, , ko, ki, . . . , k E 
{1, . . . , n}, and < < 1, /c G {1, . . . , ?7,}\{fco, ki, . . . , ki}. Without 
loss of generahty we may assume that ko < ki < ■ ■ ■ < ki. Then the 
Markov chain {x2{Tm))m>i has a unique communicating positive- 
recurrent class {(fko, fko+iy • • • ! V^fei}- All other states are transient. 
So there exists a unique invariant distribution (vr^)^^]^. Besides, it 
is not difficult to verify that 

TTfc = 0, 1 < k < ko OT ki < k < n, 
TTfco+j = fi{qko+i, ■ ■ ■ ,qki-i,Pko+i, ■ ■ ■ ,Pki-i,ki — ko), 
i = 0,l,...,ki-ko. 

where fi,i = 0,1, ki — ko, is defined by formulas (!30l) - (133|) . 

The case of pfc^ = and < pfc < 1, k E {I, . . . ,n}\ {ko} can be 
treated analogously. 

3.4. Suppose there exist ki, /c2, • • • , ^2/ ^ ■ ■ ■ , n}, ki < k2 < 
■ ■ ■ < k2i, such that p^^ = p,,^ = ■ ■ ■ = pf,^^ _^ = 1, p^^ = = . . . = 
Pk2i ~ 0' < pfc < 1, A;2j-i < k < k2j, j = I, . . . ,1. Then each 
of the classes {(fki,fki+i, ■ ■ -^fki}, {v^fcai-n V^fcai-i+i, • • • ^V^fcaJ is 
a communicating one for {x2{Tm))m>i- (Elements of different classes 
are not communicate.) For each class, there exists a unique invariant 
distribution which can be found by formulas fl30l) - fl33l) . 

Assume X2{ti) = ipk^, k2j < ko < A;2j+i for some j G {!,...,/}. 
{k2i+i '■= ki). Define 

a{ko) = P^,Jx2(a^2, A cr^2j+i) = V2j}- 

Note that 

P'/'.oi^^a, A(T<^2^.+, < 00) = 1. 
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If < pfc < 1, k2j < k < 2j + 1, computations lead to the formula 



Q!(fco) - 



1 _|_ v^fi:2j + l 



If < Pfc < 1, A;2j < ^ < ^2j+i5 and equality is reached, aiko) = 1. 
Similarly, if < < 1, A;2j < < ^2j+i, and equality is reached, 
q;(/co) = 0. 

If there exist ii, . . . , is, is+i, ■ ■ ■ such that k2j < h < ■ ■ ■ < is < 

fco < is+l < ■■■ < ir < k2j+l, Ph = ■ ■ ■ = Pis = 1, Pis+i = ■■ ■ = Pr = 

0. Then 

4. The main result. 

The behavior of the process {x{t))t>o is uniquely determined by that 
of the process {x(t))t>o. As we know (see Remark [2]) t] = oo a.s. So 
x{t) — > as t ^ (— if and only if x{t) +oo as t — +oo. 

Let us obtain an estimate of the local time L^''(t), k = 1, . . . , n. 
Fix k G {1, . . . , n}. Then the functional L^^{t) increases at those and 
only those intervals of time [r^, r^+i), > 1, for which x{Tm) € c^. 
Moreover, 

Ltit) = LZit) =■ ljm^/*2[^,„,,^,+,](x2(s))rfs. 

Let z/fc(t) be the number of such intervals on [0, t\. According to the 
strong Law of large numbers 

KkLZ((^Vk-i A (^vk+,) as t ^ +00 a.s. 

Fix e > 0. Then there exists tkiiyj) > such that for all t > tki 
the inequality 

z/fc(t)Mfc(l -6)< Lfit) < i/fc(t)Mfc(l + e) (34) 

holds. Recall that := E<^j.L|'^(cr^^_^ A c^i^^.^) and it is defined by 
formula 
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Further, by the renewal theorem (cf. [6], Ch. XI), 



t 



Kill, as t ^ oo, 



where K is some positive constant. Using the strong law of large 
numbers we get 

Mt) J. . 

> K-Kh as r ^ oo a.s. 

t 

Thus there exists tk2 > tki such that for all t > tk2 

Kuktil -e) < Ukit) < Kiiktil + e) a.s.. (35) 

Finally from (!34ll and (!35l) we obtain the inequality 

{l-efK'KktMk<L'J:{t)< 

< (l+e?K7iktMk, t>tk2, a.s. 



(36) 



We use the inequality fl5^ to get an estimation of the process 

{Xl{t))t>0. 

Let J^k^ilk'^kMkta.'nOk 7^ 0. Integrating (fT9|) over [0,t] and 
applying (!36|) we come to the inequality 

n 
k=l 

<Mt)< (37) 

n 

< xi(0) + wi{t) + KtY,lk{l+ Skf-KkMk tan Ok 

k=l 

valid for all t > t* = max^„Y^{ti2, . . . , t„2}, where 



e if 7fctan6'fc > 0, 
-e if 7fctan6'fc < 0. 



So if ^2=1 Tfc^fc^fc tan^fc 7^ we can choose e so small that the 
expressions 



^7fc7rfcMfctan Ok, 



k=l 



16 



^ 7fc(l + £^fc)^vrfcMfc tan 6'fc, 

k=l 

and 

n 

^7fc(l - £fc)^7rfcMfctan6'fc 

k=l 

have the same sign. By the low of iterated logarithm there exists 
t^: > such that for all t > t^, the inequality 

\wi{t)\ < (1 + £)V2tlnlnt (38) 

fulfilled. Combining ( 1371) and ( 1381) we see that 

— s> +00 as t 

Xi{t) — oo as t 

Now we can return from the process {x{t))t>o to {x{t))t>o making 
the inverse change of time and doing the rotation of coordinate 
system inverse to that defined by ([7j). Then we see that x{t) 0, 

as t ^ if ELi(Pfc - ^'^^'^'^ Pkt-i+qA ^^"^^k > 0. In this case, 
as a consequence of (1571) . there exist tQ,Ki, K2 > such that for all 
t > to the inequality xi(t) > Ki + holds true a.s. This implies 
convergence of the integral e~'^^^^^^ds. Therefore, C < 00 a.s. 

Similarly, if ELi(Pfc " ^^^^'^ < 0' = 

a.s. and x{t) ^ 00 as t ^ 00 a.s. 

Let Ylk=i 'y kT^kMk tan 6k = 0. Define 

Hi = mi{t > : x{t) G Ci}, 

and, for j > 2, 

= inf{t > : G Ci}, 

where 

jlj = inf {t > fij : G {c„, C2}} 

(see Fig. [3]). 



-00 if y^7fc7rfcMfctan6'fc > 

k=l 
n 

-00 if y^7fc7rfcMfctan6'fc < 0. 



k=l 
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Figure 3. 



Denote 

n 

i^j = J]7fctan4(4'=(/i,+i) - Lfifij)), j > 1. 

k=l 

It is easily seen that if Ylt=i Ik^^k^k tan 9^ = then Kipj = 0, j > 1. 
From flTIJl) we get 

j 

Note that xi = uii{fii-^.i) — Wi{fii) + ipi, I > 1, are independent 
identically distributed random variables (i.i.d.r.v.), E^; = 0, and all 
the moments are finite. Then by law of iterated logarithm for sum 
of i.i.d.r.v. (cf. [7j) we have 

lim Xi{fij) = +00, lim Xi{fij) = — oo a.s. 

It is easily seen that in this case the integral Jp°° e~^^^^*-'(is diverges 
a.s. Therefore ( = +oo a.s. 

Recall that conventionally Cni+k = Ck, ^ni+k = ^k and put addi- 
tionally -fnl+k = Ik, Tfnl+k = TTfc, k = 1, . . . , U, I E N. 

Definition 1. We call qOcj, 1 < i < j, a. characteristic wedge 
if (fi,(pi^i, . . . ,(fj is a communicating class for the Markov chain 

{x2{'Tm))m>l- 
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Theorem 1. Let x{0) = x° 7^ 0. Then 

1) P^JC < 00} + P,.JIh^^t_+oo \x{t)\ = 00} = 1. 

2) The probability Fxg{( < 00} is equal to the probability of hitting 
some characteristic wedge Cj-^Ocj^, 1 < ji < j2, for which 

E{pk — Qk)'^k — ; 7- tan Ok > 

^—ji 

where tt^, ji < k < j2, is defined in Section 3. 

In particular, suppose < < 1 or < p/t < 1, = 1, . . . , n. 
Then, if Yl=i{Pk - gfc)^^ ^^^'-^^^^- ^ ^' ^^(^ V™^^^^^ {.x{t))t>o 

hits the origin a.s.; if Y2=i{Pk - Qk)T^k^J^f^:^^tan9k < 0, the 
process {x(t))t>o does not hit the origin a.s. and \x(t)\ 00 as 
t +00 a.s.; ifY.k=iiPk - 1k)'^k pJl;~^ll^^^ tangfc = the process 
does not hit the origin a.s. and lim(_>+oo \x{t)\ =00 a.s. 

Remark 4. It is not difficult to calculate that if < < • • • • ■ 
Pn ^ Qi ■ ■ ■ ■ ■ Qn, the process {x{t))t>o hits the origin if and only if 

n n 

X] ^fc X] '^'=+1 ■ • • • ■ ''^k+j-iirk+j - 1) tan6'fc+j > 
k=i j=i 

where Vj = Pj/qj, j = I, . . . ,n. 
5. Examples 

Example 1. Let n = 1. This is a degenerate case in which the 
process reaches the origin if and only if tan^i > 0, i.e. 61 > 0. 
Example 2. Let n = 2, ^1 = ^, < pi,p2 < 1. Then tti = tts = 1/2 
and the process hits the origin if and only if 

(Pi - Qi)iP2 - g2)^(tan6'i - tan 6^2) + 

+ 271 {{pi - gi)g2tan6'i + {p2 - g2)Pitan6'2) > 0. 

Consider the particular case when pi = 1, p2 = 0. Then the part of 
the process in the interior of the wedge Ci0c2 is a Brownian motion 
with instantaneous reflection at the boundary of the wedge (Fig. 
Hj). Let the starting point x^ be inside the wedge and away from 
the origin. In this case the process hits the origin if and only if 
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Figure 4. 



tan^i — tan ^2 > 0, i.e. 6i — 62 > 0. Thus the result we obtained 
coincides with that of [1]. 

Example 3. Let n = 3, < < 1 or < < 1, /c = 1,2,3. If 
P1P2P3 7^ (liQ2(l3, the invariant distribution has the form 

q2 + P2P3 

^2 = 

qi + P1P2 
vrs = -^—^ 

where D = qi + q2 + q3 +P1P2 +P2P3 +PiP3- In this case the process 
hits the origin if and only if 

(Pi - qi) tan 6*1(52^36 + 92^36 + ^2^36) + 

+ iP2 - q2) tan6'2(pig3^3 + ^i^sG + ^1^36) + 

+ {P3 - q3) tan6'3(gip26 + ^1926 + P1P26) > 0. 
If P1P2P3 = qiq2q3 we get 

52^3 

q2q3 + Piq3 + P1P2 
Pih 

7r2 = - - - - - - , 

g2g3 +pig3 + P1P2 
P1P2 

TTs = -r-: 

^2^3 +Piq3+PlP2 

and the process hits the origin if and only if 

q2q3{Pi - gi) tan6'i +pig3(p2 - ^2) tan6'2 +PiP2(P3 - ^3) tan 6*3 > 0. 
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